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We consider interfacial phenomena accompanying bulk quantum phase transitions in presence 
of surface fields. On general grounds we argue that the surface contribution to the system free 
energy involves a line of singularities characteristic of an interfacial phase transition, occurring 
below the bulk transition temperature T c down to T = 0. This implies the occurrence of an 
interfacial quantum critical regime extending into finite temperatures and located within the portion 
of the phase diagram where the bulk is ordered. Even in situations, where the bulk order sets in 
discontinuously at T — 0, the system's behavior at the boundary may be controlled by a divergent 
length scale if the tricritical temperature is sufficiently low. Relying on an effective interfacial model 
we compute the surface phase diagram in bulk spatial dimensionality d > 2 and extract the values 
of the exponents describing the interfacial singularities in d > 3. 
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I. INTRODUCTION 

The theory of interacting quantum systems undergoing 
bulk ground state phase transitions is presently an arena 
of fascinating developments strongly linked to parallelly 
progressing experimental work* 1 - - — Most of the theoreti- 
cal interest is concerned with bulk properties of homoge- 
neous systems. Realistic experimental samples are how- 
ever finite and inhomogeneous at the boundaries. The 
issues concerning the boundaries' influence on bulk criti- 
cal behavior and the character of interfacial singularities 
at bulk criticality were widely investigated for classical 
phase transitions^ So were interface unbinding transi- 
tions, presenting a separate class of critical phenomena 
not related directly to singularities in the bulk free en- 
ergy j»— However, the problem of interfacial quantum 
critical behavior has so far received surprisingly little at- 
tention^i - — considering its potential practical relevance 
and the fact that the critical behavior at the boundary 
is completely distinct to its bulk counterpart. In this pa- 
per we adapt an argument due to Cahn, 1 - well known 
in the theory of interfacial physics, to argue that the on- 
set of bulk order at zero temperature is accompanied by 
an interfacial phase transition, which may be of second 
order even if the bulk transition is discontinuous. This 
surface phase transition can be described by a quantum 
version of the capillary-wave Hamiltonian. We subse- 
quently apply a functional renormalization-group (RG) 
technique to analyze the interfacial singularities along the 
zero isotherm, as well as approaching the quantum phase 
transition from the finite temperature region, and char- 
acterize the surface quantum critical regime. We assume 
that the underlying microscopic interactions are short- 
ranged and the bulk phase transition can be described 
by an effective bosonic d-dimensional field theory char- 
acterized by a (bulk) dynamical exponent z. This may 
apply to systems of static spinsi as well as itinerant Fermi 
systems i 15 ' 16 We focus on the case of discrete Ising-type 
symmetry breaking in the bulk. The generic bulk phase 
diagrams of the systems in question are schematically de- 



picted in Fig. 1. The transition can be tuned by varying 
a non-thermal control parameter S, or the temperature 
T. At T — there is a quantum phase transition of sec- 
ond (A) or first (B) order. In the latter case a tricritical 
point generically occurs at T tr i > 0Ji~— 
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FIG. 1: A schematic illustration of the bulk phase diagrams. 
The transition between the symmetric (SYM) and symmetry- 
broken (S-B) phases can be tuned by varying a nonthermal 
control parameter 5 or the temperature T. The second-order 
phase transition line (T c -line) terminates at a quantum critical 
point (QCP) (A), or becomes first-order below a tricritical 
temperature Ttri (B). The first-order transition line below 
the tricritical point (TCP) is represented (schematically) by 
the vertical green line. 

In the present work we consider the system in the 
symmetry-broken phase (T < T C (S)), and subject to an 
infinitesimal bulk field h — 0~ favoring one of the two 
available bulk phases (phase a). On the other hand, we 
impose an external surface field h± > acting only at 
the boundary of the system and giving rise to an oppo- 
site ordering tendency (phase /?) at the surface. This 
implies the formation of a layer of the /3-like phase at the 
wall. The layer may be either microscopic or macroscopic 
in width, the two situations being separated by what is 
called an interface unbinding or wetting transition ) 14 : 23 
at which the a — f3 interface becomes macroscopically sep- 
arated from the wall, and the capillary fluctuations di- 
verge. Such a phenomenological picture may be realized 
in microscopic models as well as at the level of coarse- 
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grained Landau-type theories supplemented with surface 
terms X Another approach employs effective interfacial 
models, which derive from further coarse-grained Lau- 
dau theoryi^ir— The classical versions of all these were 
extensively studied in course of the last 35 years. The 
quantum variant, relevant for the case of temperatures 
approaching zero seems however missing. 
Particular questions addressed in the present paper are: 
(1) Is the shape of the quantum interface unbinding tran- 
sition line universal in the limit T — > like happens to 
be in the bulk casej 6 ' 16 i 27 (2) Is there an analog of the 
bulk quantum critical regime associated with the quan- 
tum interface unbinding transition, if so, what are the 
properties of the system in this regime; (3) The classi- 
cal interface unbinding transition in d = 3 is extremely 
unusual^ - — being characterized by nonuniversal crit- 
ical exponents and in certain cases essentially singular 
behavior. Does this exotic, nonuniversal behavior persist 
at T — > 0, how is this altered by quantum fluctuations? 

The paper is structured as follows: In Sec. II we ap- 
ply a standard argument due to Cahn to argue, that bulk 
quantum criticality occurring in presence of suitably cho- 
sen surface fields is generically accompanied by a separate 
phase transition restricted to the system boundary. In 
Sec. Ill we introduce a phenomenological effective Hamil- 
tonian model capturing the relevant thermal and quan- 
tum interfacial fluctuations. In Sec. IV we derive the 
functional RG flow equation to be solved in Sec. V at zero 
temperature and in Sec. VI in the classical limit (stud- 
ied already previously). Relying on the results of Sec. V 
and Sec. VI we are able to approximately solve the RG 
flow equation at finite but low temperatures, traversing 
both the T = and the classical scaling regimes. This is 
achieved in Sec. VII yielding results summarized in phase 
diagrams. We conclude in Sec. VIII. 

II. CAHN ARGUMENT 

The original argument due to Cahn— refers to a situ- 
ation, where a first order classical phase transition line 
terminates at a classical critical point. The temperature 
is increased so that the system approaches the critical 
point, remaining infinitesimally close to the transition in 
the phase disfavored by the boundary. An interface un- 
binding transition must occur at the wall before the bulk 
critical temperature is reached. We present the argument 
in a version adapted to the presently relevant setup. We 
approach a second-order bulk phase transition (Fig. 1) 
along any isotherm T > (case (A)) or T > T tr i (case 
(B)). The bulk field h — CP and the surface field is fixed 
at a positive value h\ > 0. Assuming no macroscopi- 
cally wide layer of the /3 phase intrudes the vicinity of 
the wall, the surface free energy density a wa and the free 
energies (per unit area) associated with the wall-/? and 
a-/3 interfaces (cr w p and a a p respectively) must fulfill the 
relation 

o wa — a w /3 — a a /3 < , (1) 



otherwise the system would find it favorable to unbind 
the a — f3 interface from the wall and form a microscop- 
ically wide layer of the (3 phase. In the vicinity of the 
critical line however cr Q( g ~ (<5 C — 8)^ and \cr wa — a w p\ ~ 
(S c — S)^ 1 . Here fli is the exponent describing vanish- 
ing of the order parameter at the wall as the critical 
point is approached^ The function 8 C (T) describes the 
line of critical points and T c (8) is its inverse. Plugging 
in the values relevant for the d = 3 Ising universality 
class (/j, s» 1.3, j3\ ps 0.8) j& one realizes that condition 
dU) is violated sufficiently close to the critical line T c (8), 
implying that at S W (T) < S C (T) the inequality ([T]) is re- 
placed with an equality and the interface unbinds. Here 
S W (T) (or T w (5)) denotes the interface unbinding transi- 
tion line. The conclusion also holds if mean-field values 
of the exponents fi and fix are used, as well as for the case 
T = T tri l For this reason the T^-line cannot end at the 
TCP. For a recent study of classical wetting in systems 
involving a tricritical point see Ref. l3li The argument 
applies also for T = in case (A), if the bulk transition 
at T = falls into the mean-field universality class, which 
is the usual situation and which we here assume. This 
also naturally comes out in the calculation in Sec. V. As 
a consequence, a layer of the j5 phase must be coating 
the wall for 8 e [8 W (T), 8 C (T)[. In case (B) the argument 
can be repeated approaching the trictitical point along 
the triple line and two situations must be distinguished 
depending on whether the interface is unbound down to 
T = or not (see Fig. 2). The former offers an interest- 
ing route of realizing quantum critical behavior at the 
boundary in the absence of any bulk quantum critical 
point. As is well known, the interface unbinding tran- 




FIG. 2: A schematic illustration of the phase diagrams as 
implied by Cahn's argument. A line of interface unbinding 
transitions necessarily continues to T = if the bulk phase 
diagram features a quantum critical point (A). In the case of 
a first-order bulk phase transition below T tr i two possibilities 
occur: The T^-line continues to T — (Bl), or terminates at 
the finite-T bulk phase transition line (B2). 

sition is accompanied by a divergent correlation length 
£|| controlling the capillary-wave like fluctuations of the 
a — ft interfaced On the other hand the bulk correlation 
length £{, remains microscopic at 8 W (T). 

Let us also note that the Cahn argument breaks down 
in a number of physically interesting cases involving long- 
ranged microscopic interactions. 32 However the quantum 
systems we address here either interact via short-ranged 
forces (static spins) or display efficient screening mech- 
anisms (itinerant fermions). We expect that the Cahn 
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argument is generically valid in the present context. 

III. MODEL 

In this section we introduce a phenomenological 
capillary-wave type model adapted to a situation where 
quantum-mechanical and thermal fluctuations must be 
treated on equal footing. The model retains the posi- 
tion and imaginary-time dependent location Z(x, t) of the 
a — (3 interface above the wall as the fluctuating degree 
of freedom. It relies therefore on the assumption that all 
the soft modes at the interface unbinding transition are of 
capillary-wave type or that the effect of other soft modes 
amounts to renormalization of the parameters describ- 
ing the effective order-parameter action. In the classical 
case this may be proven in d = 2£2- The quantum ef- 
fects are incorporated into the model via the imaginary 
time dependence of the fluctuating degree of freedom. 34 
In the classical context the capillary-wave Hamiltonian 
was derived 2 ^ starting from an underlying bulk Landau 
theory supplemented by surface terms. In such derivation 
one fixes the position of the phase interface, constraining 
the accessible order-parameter profiles, and subsequently 
performs a partial trace respecting this constraint. Since 
no soft modes responsible for the critical singularities 
are integrated out in this procedure, one may rely on 
a rather crude approximation retaining only the domi- 
nant (saddle-point) term in the constrained trace. For 
the present case, involving quantum-mechanical effects, 
such analysis may be carried out along fully analogous 
lines. For systems of static spins the relevant effective 
bulk theory is a <i+l-dimensional <f> model, where finite- 
temperature effects are accounted for by a finite-size ef- 
fect in one direction^ The analysis of Ref. [26| general- 
izes to the quantum case straightforwardly and yields 
the interfacial critical exponent z = 1. The situation is 
somewhat more complicated in Fermi systems, charac- 
terized by a Fermi surface and soft bulk excitations oc- 
curring in the entire thermodynamic phases. This yields 
the bulk dynamical exponent z > I. 15 For such cases, 
in the present study, we shall assume that the damping 
term in the effective action is inherited from the bulk 
order-parameter theory. We shall therefore neglect the 
possibility that direct coupling of the fluctuating inter- 
face to the soft particle-hole excitations across the Fermi 
surface changes the way in which the imaginary time oc- 
curs in the effective action. For such cases we shall focus 
on z = 2 pertinent to itinerant Fermi systems ordering 
at wavevector Q ^ 0, which may apply either to charge- 
or spin-density-waves. 

Relying on the above we consider the bare effective 
action: 

W[I(x,t)] = j ' dtj d d - x xV{l) 

+T E/ (l^r^ 2 + ^M 2/ 1W- q - • ( 2 ) 



Here V(l) is a local potential describing interactions be- 
tween the wall and the a — /3 interface, a is a constant 
interfacial stiffness coefficient, Z qtJ is the Fourier trans- 
form of Z(x,t), and z € {1,2}. Finally w — 2imT 
are the (bosonic) Matsubara frequencies, and the quan- 
tum damping term involves a constant Z u . As already 
stated, this model may be understood as a coarse-grained 
Landau-Ginzburg Hamiltonian, out of which the bulk or- 
der parameter fluctuations had been integrated out, re- 
taining £(x, t) as the only degree of freedom. The inter- 
facial potential reads 

V(l) =re-^ +be- 2l/ ^ + ... . (3) 

The parameters r ~ (T — T w ) and b > may (alike in 
the classical case) be expressed in terms of the underlying 
Landau-Ginzburg model. For 2 = 1 pertinent to systems 
of static spins the model defined by Eq. (2) and Eq. (3) 
may be derived from the bulk order-parameter action 1 
parallelizing Ref. [26|. Indeed, the quantum effective ac- 
tion is in this case equivalent to the standard classical 
<^ 4 -model in dimensionality elevated by 1, and where the 
system size in one direction is p = 1/T. For this rea- 
son it is also very plausible that the capillary waves are 
the only relevant low-energy modes also for the quantum 
systems. On the other hand, the value of z is usually 
altered in the cases where the order parameter field in- 
teracts with a bath of generically soft modes like happens 
to be in the Fermi systems. In such cases one may expect 
that the value of z depends on the particular bulk insta- 
bility in question and potentially other characteristics of 
the system. We have performed explicit calculations for 
z = 1,2. As the calculation below shows, in d > 3 the 
interfacial quantum critical behavior does not show a de- 
pendence on the value of z apart from coefficients which 
are nonuniversal anyway. 

We skip higher-order terms and the dependence of a on 
in Eq. ([3]) At mean-field level the average separation 
of the interface l„ from the wall results from minimizing 
V(l). The interface unbinds as r ~ (6 — S W (T)) -» 0~. 
One straightforwardly obtains l n = log and the 
interfacial correlation length £jj~ 2 ~ V"^) = r 2 /(26). 
Note that for t > the interface remains unbound 
and the system features soft capillary fluctuations with 
= 0. The quantities r and b > can be re- 
lated to the parameters of a "semi- microscopic" Landau- 
Ginzburg theory. In particular one finds 26 r ~ h\ + grriQ, 
where g is a surface enhancement parameter and too the 
average bulk order parameter. A much more complete 
description of classical wetting is achieved within nonlo- 
cal interfacial models^^ where the leading interfacial 
singularities however turn out to be essentially of the 
same type as those resulting from the much simpler local 
approach. In Eq. Q the variables should be understood 
as coarse-grained over a scale of the order of the intrinsic 
width of the interface, which in turn is typically of the 
order of a few bulk correlation lengths. The Hamiltonian 
Eq. @ is therefore supplemented by an upper momen- 
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turn cutoff Ao somewhat lower than £ b 1 . In the practical 
calculations we will put A = t;^ 1 = 1. 

IV. RG FLOW EQUATIONS 

Description of the classical interface unbinding transi- 
tions is among the early successful applications of func- 
tional renormalization group. 37 Within a linearized treat- 
ment the critical singularities may be captured in d > 3 
by an analytical description. 2 §±2& The present paper may 
be viewed as an extension of the studies of interface un- 
binding transitions to the case involving the limit r^O 
and requiring an account of quantum fluctuations. In 
what follows we derive the relevant RG flow equations 
from an exact functional RG flow equation for the effec- 
tive average action Ta^], which is the generating func- 
tional for one-particle irreducible vertex functions in the 
presence of an infrared cutoff at scale A. The cutoff pa- 
rameter A £ [0, Ao] and the effective action interpolates 
between the bare action (Hamiltonian) at A —¥ Ao, and 
the full effective action (Gibbs free energy) as the cutoff 
is removed A — > 0. The exact flow equation 3 -^ reads: 



1 



d A T A [l] = 2 Tl ^ (2 ) 



d A R A (q) 



TX'[l] + R A (q) 



(4) 



where R A (q) is a momentum cutoff function added to 
the inverse propagator to regularize the infrared diver- 



gencies, and T^f/] 



The trace 



sums over momenta and Matsubara frequencies: Tr = 
T I (t^- 1 ■ Integrating Eq. (U]) with the initial con- 
dition at Ao specified by Eq. ([2]) would yield full solution 
to the model. Herein we will approximate T A [l] by the 
following ansatz: 



l/T 



dt I d d - 1 xV A (l) 



-T 



Jd-l, 



T^tN 2 + Z^\uj\ 2 I z \1^1 



q— uj 



(5) 



This approximation retains functional flow of the local 
potential V\(l) and completely neglects the renormaliza- 
tion of the propagator. The latter does not influence 
the critical singularities, since the anomalous dimension 
rj for interface unbinding transitions is robustly zero 3 ^ 
and, as we show later on, the precise value of z is of mi- 
nor importance to the critical behavior. The fact that 
rj = Q makes the present approach easier than an analo- 
gous treatment of bulk criticality, where anomalous scal- 
ing of the propagator requires account of its flow under 
the RG transformation.— ~— Evaluation of Eq. (J4j) for a 
uniform field I yields the flow of V(l): 

q T/ m l rr d A R A (q) , . 

where primes denote differentiation with respect to I. In 
the classical limit, where only the contribution from lo — 



is retained, Eq. (J6)) is scale invariant when rewritten in 
the variables 

A q 2 ~ _ 3^ 

s = - log j- , y = > b = v a A 2 I , 

v A (b) = V A (l)A- d+ \ r(„) = ^. (7) 



aq" 



Namely, we obtain: 



Q 7 

d s v A (b) = (d- l)tt A (6) + -^-bd- b v A (b) 



S, 



d-2 



4(2tt 



id— 1 



dyy~ 



2 

d s r A (q) 



1 + y^dfv^b) + r A (y) 



(8) 



where we performed the angular integration, and Sd — 
2 7r (rf+ 1 )/ 2 /r((d + l)/2) is the surface area of a d- 
dimensional sphere of radius 1 . In the sharp cutoff limit 



R A (q) = lim ctAV(A 2 - 9 ) 



(9) 



this is equivalent to the continuous limit of the Wilson's 
approximate recursion relations as derived and analyzed 
for classical interface unbinding phenomena in Ref. 130 . 

On the other hand, in the zero-temperature limit 
T — » we have TV„ -> f ^ and the scale invariance 
1 : Eq. ([SJ is achieved by the transformation Eq. ([7} after 
substituting d — > d + z. This is a well-known fact, that 
many quantum-critical systems (at T = 0) fall into the 
same universality class as their classical counterparts in 
dimensionality elevated by the dynamical exponent. This 
is also true for the interfacial phase transitions considered 
here. 

In what follows, we shall rely on a linearized version of 
Eq. ([6]) . Such an approach exploits the fact that close to 
the transition the effective potential V(l) is very fiat ex- 
cept for the vicinity of the wall, and the minimum of V A (l) 
is well separated from the system boundary {l^/^b 3> !)■ 
The linearized treatment captures the leading classical 
critical behavior in d > 3. After expanding Eq. (|6]) and 
dropping a constant term which merely shifts the effec- 
tive potential without influencing its shape, we obtain 



1. 



d A R A (q) 



12™ 



Z^l* + R A (q)Y 

(10) 

We note, that in the classical limit, for d = 3 and after 
rescaling according to Eq. (J7J this equation is equivalent 
to the linear RG flow equation analyzed in Ref. [29|. This 
result does not depend on the particular choice of the 
cutoff function R A (q)£l We rewrite Eq. (fTUj) as 



d A V A (l) = -/ T (AK'(Z) , 



(11) 



and split the RHS /t(A) = f#(A) + /|(A) into a clas- 
sical contribution /j?(A) from uj = 0, and a quantum 
contribution /|>(A) from lo ^ 0. 

Following, we exploit the freedom of the cutoff choice 
by taking^ 



R A (q) = a(A 2 - q 2 )9(A 2 - q 2 ) 



(12) 
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With this choice the q-integration in Eq. (fTTj) becomes 
restricted to q < A and aq 2 + R\(q) — <jA 2 for q < A. 
We obtain 

cd— 2 TiAd-4 

^ = Wy^T^W, (13) 

The Matsubara sums in f T (A) are preformed using stan- 
dard techniques. 
For z — 1 we find 



s- 



d-2 



TaA a 



ir 4 f' 4 x 



f ^ (2n) d - 1 2{d - f )Z2(2ttT) 4 ' 
-2 + T"- 1 coth(T"- 1 ) +f'- 2 sinh" 2 (f'- 1 )) , (14) 



while for z 



/1(A) - 



S d ~ 2 To-A d 2*i(l+f- 1 ) 
(2 7 r) d - 1 rf-l (27rTZ u ) 2 ' 



(15) 



where the polygamma functions ^k(z) are defined via 
* fc +i(«) = fc = 0, 1, 2, ... and = T'(z)/r(z). 

We introduced f = and f" = For future 

reference we also define T z — 2ir z ~ 1 (% i ) Z ^ 2 tj- 

Since we are interested in the low temperature regime, 
T z is small at the beginning of the flow, and the term in- 
volving f T (A) dominates. As the cutoff scale is reduced, 
T z grows and the flow may cross over to the classical 
regime dominated by the term with /j? (A). This will 
however clearly never happen at T = 0. 

Eq. (|11|) can be integrated starting from an initial po- 
tential Va o (0- This yields 



V A (Z) = I ^'(-4 7 r J F T (A))- 1 / 2 e*™^A a') . (16) 



where 



F T (A) 



Ao 



dA'MA') 



(17) 



We note that although the flow equation was derived 
from an exact functional flow equation, the problem at 
hand is absolutely perturbative (at least in d > 3). In- 
deed, the same linear flow equation can be derived via 
a somewhat more cumbersome traditional procedure of 
Wilsonian momentum shell integration. We also note, 
that the present approximate linear RG approach does 
not yield any fixed point except for the non-interactingi^S, 
The procedure to be applied in Sec. VII in order to com- 
pute the phase diagrams combines those of Ref. [H and 
Ref. [29]. We renormalize to a scale at which the cur- 
vature of the renormalized, rescaled potential becomes 
of the order unity (in units where the microscopic scale 



An 



1). We then treat the renormalized potential 



with a saddle-point approximation and rescale back to 
the physical lengths. The flow is split into a quantum 
(T z < 1) and classical (T z > 1) stages. We scale the vari- 
ables according to Eq. in the classical regime, and 
replace d with d + z in Eq. ([7]) for the flow in the quan- 
tum regime. 



V. SOLUTION AT ZERO TEMPERATURE 

We present details of the solution for z = 2. The 
procedure for z = 1 is fully analogous. At T = 
the classical term f T l (A) — and T <C 1 at any A. 
By expanding the polygamma function around T = 0: 
#i(l + f- 1 ) = f + 0{f 2 ) we find 



/ 9 (A)=^im/!(A) = C d A 



d-2 



(18) 



where we introduced C a 



25" 



The integral 



^(d-l)(2ir) d ' 

defining Ft (A) is now straightforward and one finds 
V A (l) = / dl'(AnK(A))-V 2 e-V- l 'f/* K Wv Ao (l>) , 



(19) 

where K(A) = - A d_1 ). This implies that the ef- 

fect of fluctuations at any d > 1 is merely to smear the 
bare potential by convoluting it with a Gaussian of finite 
width. This is reminiscent of the solution to the classi- 
cal problem in d > 3^ and consistent with the expec- 
tation that the quantum interface unbinding transition 
at T — is in the same universality class as its clas- 
sical counterpart in higher dimensionality D = d + z. 1 
Note that the transition occurs at r = and the crit- 
ical value of S W (T — 0) is not affected by fluctuations. 
In the limit t -> 0~ (S -» S W (T = 0)~) one finds 
^ij" 1 ~ (5 W (T = 0) — S) and a logarithmic divergence 
of l n . 

The solution proceeds along the same line for z = 1, 
where in Eq. (fT5|) A d ~ 2 is replaced by A d ~ 3 and the con- 
stant C d takes a different form. Again, the renormalized 
potential V\(l) is nothing more than a convolution of the 
bare one V Ao (I) with a Gaussian whose width is of micro- 
scopic extent and the critical properties are described by 
the mean field theory (MFT). The general solution for 
z G {1,2} can be written as 



V A (l)= I di'-J—e-C'-WAiv, (0 , 

V nA A 



(20) 



where 




{2ir) d Z ul (d-l) 



VI. SOLUTION IN THE CLASSICAL LIMIT 

As we already mentioned, in the classical limit the 
quantum contribution to the flow equation vanishes 
(f T (A) — > 0), the flow is dominated by the classical term, 
and up to a multiplicative constant A(d) the linearized 
flow equation is equivalent to the equation analyzed in 
Ref. [29|. Occurrence of A(d) is due to the cutoff choice, 
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which is different here as compared to Ref. [29|. Impor- 
tantly A(3) = lr^ which holds irrespective of R A (q). On 
the other hand, above the upper critical dimension d > 3, 
the quantity A(d) has no impact on the critical behav- 
ior. Our results in the high temperature limit coincide 
with Ref. H^, and we quote them here. The bare effective 
potential 



V Ao (l)=A Ao (l) + R Ao (l) + W Ao (l) 



(22) 



is split into an attractive tail A Ao (l) = Te~ l '^ b 6(l), a 
repulsive tail R Ao (l) = be~ 2l ^ b 9(l) and a wall repulsion 
contribution W Ao (l). Following Ref. 1 23 we model the lat- 
ter by a soft- wall potential W Ao (l) = w8(—l), w being a 
constant. For d > 3 fluctuations do not introduce sig- 
nificant renormalization to the bare potential, and the 
renormalization amounts to convoluting the bare poten- 
tial with a Gaussian of finite width. The critical singular- 
ities are identical to those of MFT. The situation is quite 
distinct in d = 3, where the convolution width becomes 
infinite. In particular, for a range of parameters, the wall 
repulsion renormalization becomes essential for the criti- 
cal behavior. The three contributions to the renormalized 
potential from the different interaction terms in Eq. (|2"2")l 
must be carefully analyzed at the scale A* , where the cur- 
vature at the minimum is of the order unity. Applying 
a mean-field approximation at this scale and returning 
to the physical length units yields the critical behavior. 
The interplay between the three terms in Eq. (|22p leads 
to three fluctuation regimes showing non-universal criti- 
cal behavior determined by the dimensionless parameter 



T 

47T<7 



(23) 



For the capillary correlation length exponent v» defined 
by £|| ~ (S w — 5)~ u w one finds 



(1 - w™)- 1 



for lj w G (0, 1/2) 
- 2 for u w G (1/2, 2) 
for oj w G (2, oo). 



For uj w — > 2 the exponent v\\ diverges and one finds 
an essential singularity, whereby £m ~ e a /i 5 ^~ s )_ The 
film thickness diverges logarithmically for lj w < 2 and 
1^ ~ 1/(S W — 5) for lo w > 2. Remarkably, the critical 
value S w is not shifted by fluctuations for u w < 2. One 
of the major questions raised in this paper concerns the 
fate of this extremely unusual picture in the quantum 
limit. 



VII. COMPUTATION OF THE PHASE 
DIAGRAM 



We now present a calculation of the interfacial phase 
diagram at T > 0. We integrate the linearized RG flow 
Eq. (TTTj) at finite, low temperature (such that T z <C 1 at 



A = Ao). Eq. (|22| acts as the initial condition. The ini- 
tial stage of the flow is dominated by the quantum contri- 
bution and we expand Eq. (fTTj) around the T — > limit, 
and rescale according to the zero-temperature dimen- 
sions, where r = A u A {r) = V A (l)hr d ~ z+1 . For 
simplicity in the definition of r we dropped the \fo factor, 
which in the present context is just a scale-independent 
constant anyway. Following the flow of the renormal- 
ized, rescaled potential, we identify two regimes in the 
space spanned by (S,T). In the first case, the parameter 
T z < 1 at the scale A*, where the potential curvature 
at the minimum becomes of order unity. We terminate 
the flow at A* and, since the scale where (A) becomes 
significant was never reached, conclude that the system 
properties are essentially the same as in the T = case. 
This regime is an analog to the quantum (or T = 0) 
regime in the case of the well-developed theory of bulk 
quantum phase transitions. The other regime occurs 
if the scale A* is not reached before T z m 1, whereby 
A = A cr . At A < A cr the flow becomes dominated by 
the classical term (A). We now neglect the quantum 
contribution /|,(A), and perform the classical rescaling. 
The flow equations become equivalent to those occurring 
in the high-temperature limit, where A c ; acts as the up- 
per cutoff and the initial condition is provided by V Acr (I) 
obtained by the first stage of the flow integration. The 
(approximate) matching procedure at A cr is analogous 
to that of Ref. [HI, and the regime A < A cr is studied 
along the same lines as Ref. [29|. We give details of the 
calculation below, and the results are summarized in the 
phase diagram Fig. [3l 



A. Flow for A > A cr 

According to the procedure outlined above, we calcu- 
late the contributions to the renormalized potential V A (l) 



from the three terms in Eq. (f2"2l . We put A 
and consider r < 0. We obtain: 



A A (l) = I re VAl7(I^)-* Erfc 



W A (l) = -wErfc 



A1 



21 



= 1 



(24) 



(25) 



(26) 



We observe that remains small and only the terms 
involving I in the arguments of the error function con- 
tribute to the renormalized potential at large I. The term 
W A (l) decays as e~ l at I 3> 1, and rapidly approaches w 
for negative I. Expanding the error functions yields: 



A A (l) 



re 



y/Af/T&r-l 



9(1) 



(27) 
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Tw-line 



Tc-line (bulk) 




QCR 
\ (bulk) 



IOCP 



OCP 



FIG. 3: A schematic illustration of the computed interfacial 
phase diagram for case (A) (see Fig. 2). The interfacial part 
of the diagram is the same for the case (Bl). The interface 
unbinding transition line T w is parallel to the bulk T c -line and 
terminates at T — with an interfacial quantum critical point 
(IQCP). The power law describing the shape of the T w -\ine 
is the same as the for the T c -line and is determined by the 
bulk properties (dimensionality and the bulk dynamical ex- 
ponent). The interfacial correlation length £m is infinite in the 
region between the T w and the T c lines. A crossover line sepa- 
rates the interfacial quantum regime (IQR) from the quantum 
critical regime (IQCR). The interfacial dynamical exponent z 
determines the shape of this crossover line, which is straight 
for z = 1 and parabolic for z = 2. In both regimes and for all 
d > 3 the interfacial height l n diverges logarithmically upon 
approaching the IQCP. The behavior of the correlation length 
£ii is always mean-field in d > 3. In d = 3 it is mean-field in 
the IQR regime, and £y ~ T lf( - x ~ DsT) when the IQCP is ap- 
proached along a straight line from the IQCR regime. 



R A (l) « beV^-^eil) 



W A (l) « w0(-l) 



(28) 



(29) 



We neglected irrelevant Gaussian tails. 

We now rescale according to r — A 2 \ y v\(r) = 
VA(Z)A~ d ~ z+1 and analyze the minimum r n of the 
rescaled flowing potential v&(r), and the curvature 
v'j^(r„). We consider two scales: A cr , where the param- 
eter T z — 1 and the flow crosses over to the classical 
scaling regime; and A*, where v'J^(r„) = 1. If A* > A cr 
classical scaling is never reached and the critical behavior 
is the same as for T = 0. Solving the inequality 



A* > A r 



(30) 



we find 



z/2 



T < (-tW- 1 (^Y b- z / 2 e z V^/(^) . (31 ) 



The central observation here is that the crossover line 
approaches r = linearly for z = 1 and parabolically 
for z — 2. The result given by Eq. (|3Tj) is valid in any 
dimensionality d > 2. Apart from the constant encoded 



in Cf the shape of the crossover line does not depend 
on the dimensionality d (which is also the case of the 
corresponding bulk system in d + z > 4) . 

For the opposite case, where Eq. (1511) does not hold, 
the flow reaches the scale A cr and crosses over to the 
classical scaling regime. 

B. Flow for A < A cr 

In the classical scaling regime A < A cr the flow be- 
comes dominated by the classical term /j? (A) in Eq. (jlljl . 
The quantum-mechanical effects are encoded in the po- 
tential V\ cr (l), which was obtained by integrating the 
quantum regime of the flow, and which now serves as 
the initial condition. Therefore both quantum and ther- 
mal fluctuations may now influence the system proper- 
ties. We consider the solution Eq. (Tl6l) to the linear flow 
equation Eq. treating A cr as the upper scale and 
VActCO as computed above as the initial condition. We 
neglect the term / T (A), rescale according to Eq. J7]), and 
terminate the flow at the scale A* , when the curvature of 
the rescaled, renormalized potential at the minimum be- 
comes of order unity. Since the correlation length at this 
scale is one, the physical correlation length £y is of order 
(A*)~i. ex t rac t information on the critical behavior 
of the interfacial height analogously. 
From the technical point of view this part of the analy- 
sis parallels Ref. [2^. The differences amount to treating 
the temperature-dependent scale A cr as the upper cutoff 
and the potential computed in the preceding section and 
evaluated at A cr as the initial condition. For the present 
case 



f 2d? (A 

-F T (A) = > rf - 3 K 



A 



d-3\ 



D d T log 



for d > 3 
for d = 3 



with 



S 



d-2 



(2 7 r) rf - 1 (d- 1)<7 



(32) 



The convolution width as given by —Ft {A) is obviously 
finite for any T > for d > 3 and the resulting critical 
behavior is purely mean-field. The situation is somewhat 
more interesting for d = 3. Parallelizing Ref. [5^ we find 
that at low T the system will find itself in a regime anal- 
ogous to the classical ui w < 1/2 case, where oj w ~ T. 
This happens because a temperature factor suppresses a 
quantity analogous to the capillary parameter u w in the 
present calculation. We find: 



log A* =log 



7>2 



l/(2-2D 3 T) 



2 J P 3 riogA c , 
2 - 2D-iT 



(33) 



where r 1 = Te V A I/ 167r , b' = beVH/*. R follows that 
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1/(1-D 3 T) 



VIII. SUMMARY 



(l- ^CfJ(M^)2^- 1 (Zjay' 2 T) 1/{1 ° 3T) . (34) 
The equilibrium interfacial width is given by 

/ ff = -logA*(l + 2D 3 T) . (35) 

Approaching the interfacial quantum critical point along 
a straight line we put r ~ — jT, 7 > 0. The conclusions 
of the calculation in d = 3 arc summarized as follows: 

1) The leading singular behavior of the correlation length 
is of the type £11 ~ j'-Vli-AsT) i m piy m g that the tem- 
perature dependence appears in the exponent. As T — > 0, 
the exponent approaches the MF value —1. 

2) The value of the dynamical exponent does not influ- 
ence either the leading singularity or the leading correc- 
tion (apart from a coefficient). This will not be true for 
d < 3. 

3) The divergence of the interfacial width in the quantum 
limit T — > is always logarithmic. Temperature depen- 
dencies due to quantum effects merely renormalize the 
proportionality coefficient. 

The calculation shows that the critical behavior occur- 
ring in d — 3 when we approach the interfacial quantum 
critical point from the interfacial quantum critical regime 
is equivalent to the case of classical transition at T 0. 
This situation is unlike the corresponding bulk transition 
and is hard to predict without performing the calculation. 
Again, we do not expect this equivalence to persist for 
d < 3. 

Neither the model defined by Eq. (2) and Eq. (3) nor 
the above RG analysis require that the bulk transition is 
second order at T = 0. Therefore our results are valid 
also for the case Bl in Fig. 2. Case B2 displays no inter- 
facial QCP and is not covered by the present caluclation. 

C. Interfacial phase transition line 

As we have observed, the interfacial transition tem- 
perature is not renormalized by fluctuations in any of 
the cases studied here. The shape of the T w -line is 
therefore fully determined by the intrinsic dependence 
of the parameter r on temperature. Invoking Ref. [26| 
we find that r ~ hi + grriQ. When the temperature 
is varied, the bulk order parameter ttlq ~ (T c — T) 13 . 
This implies r ~ (6 - 5 W {T)) = d + C 2 (S - 8 c {T)f . 
For 5 = S W (T) we find S C (T) - S W (T) = const. For 
d > 2 one has T c ~ (5 — 5q)^, where the shift exponent 
ijj = z/(d + z - 2)J£ The shape of the 7^-line follows the 
same power law as the T c -line in d > 3. 

We have however nothing to say about the shape of 
the T^-line in cases Bl and B2 (Fig. 2). Alike the first- 
order bulk transition lines, these curves are not expected 
to display universal behavior. 



We have analyzed interfacial phenomena occurring in 
presence of bulk quantum phase transitions. Applying 
a reasoning adapted from the theory of classical inter- 
facial phase transitions we have argued that interfacial 
unbinding transitions generically accompany quantum 
critical phenomena. Notably, the boundary transitions 
may occur (and be of second order) even if the bulk 
transitions are discontinuous. An interface unbinding 
transition implies a divergent interfacial length scale in 
a portion of the phase diagram within the bulk-ordered 
regime. This in turn gives rise to two scaling regimes 
in a region where the bulk is ordered but the interface 
is bound. We have found that, unlike the case of bulk 
transitions, the interfacial transition line in the (5, T) 
plane is not shifted by fluctuations, and (for the case of 
a continuous transition in the bulk) is determined by 
the shape of the bulk T c -line. Relying on a quantum 
version of an effective interfacial Hamiltonian, we have 
computed the crossover between the quantum (T = 0) 
and the quantum-critical interfacial regimes. Our result 
for the shape of the crossover line should be valid in 
arbitrary dimensionality d > 2. The critical behavior 
in the quantum regime is purely mean- field for d > 2, 
as anticipated from quantum to classical mapping 
arguments. So is the system's behavior in the interfacial 
quantum-critical regime for d > 3. Our calculations 
in the quantum critical regime are reliable only for 
d > 3. In the d = 3 case the interfacial correlation 
length diverges with a temperature dependent criti- 
cal exponent, where the temperature-dependent part 
asymptotically vanishes. This feature is closely related 
to the nonuniversality of the classical critical wetting in 
d = 3. The crossover between the two interfacial scaling 
regimes is rather weak and in fact detectable only via the 
temperature-dependence of the interfacial correlation 
length exponent mi . In addition to the interfacial cor- 
relation length, we have computed the critical behavior 
of the average interfacial width, and found logarithmic 
singularities in both scaling regimes. We anticipate that 
the crossover between the two regimes is much more 
pronounced for d < 3. In particular, for d = 2 and z = 1 
the system at T — may be in the same universality 
class as classical wetting in d = 3 yielding a very reach 
spectrum of possible critical behavior. On the other 
hand, the classical interface unbinding transition in 
d = 2 is governed by non-mean-field universal exponents. 
The properties of the quantum-critical regime in d = 2 is 
an interesting problem, beyond the reach of the analyt- 
ical linearized functional RG methodology, but open to 
be addressed within a numerical approach in near future. 
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